Perfect memristor: non- volatility, switching and negative differential resistance 
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We propose a simple model of a nanoswitch as a memory resistor. The resistance of the nanoswitch 
is determined by electron tunneling through a nanoparticle diffusing around one or more potential 
minima located between the electrodes in the presence of Joule's heat dissipation. In the case of 
a single potential minimum, we observe hysteresis of the resistance at finite applied currents and a 
negative differential resistance. For two (or more) minima the switching mechanism is non-volatile, 
meaning that the memristor can switch to a resistive state of choice and stay there. Moreover, the 
noise spectra of the switch exhibit 1//^ -^ I// crossover, in agreement with recent experimental 
results. 
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Fast progress in flash memory suggests that standard 
magnetic hard drives will soon be replaced by faster and 
much more reliable solid state drives [l[. However, solid 
state drives are still very expensive and their capacitance 
is very limited when compared with magnetic drives. To 
overcome these limitations, a new generation of flash-like 
memory has to be developed [21 • 

A most promising candidate for very fast and high den- 
sity memory are nanomemristors (memory- resistors) |3|. 
In particular, resistive switching was observed in binary 
oxides, like Ti02, NiO and perovskite oxide memristors. 
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FIG. I: Schematic diagram of the proposed memristor driven 
by either a current or voltage source: (a) A nanoparticle 
(solid red circle) can fluctuate between two electrodes (blue 
squares), due to thermal noise. This results in a fluctuat- 
ing resistance of electrons tunneling through the system (see 
text). (b),(c) The nanoparticle is subject to either a single 
well or a double well potential, mirror symmetric in (b) and 
asymmetric in (c). 



Memristive switches would provide a viable solution for 
dense memory, logic, and neurocomputing, if they can be 
demonstrated to advantageously scale with size, power, 
and driving voltage, and perform without failure and fa- 
tigue in a repeatable and uniform manner. Some of these 
issues are not well understood although the hysteretic be- 
havior of such materials, especially thin films of transition 
metal oxides like Ta205, NbzOs, Ti02, NiO, CU2O in 
metal-insulator-metal vertical devices has been known for 
decades [j|. In recent years, the interest in various oxide- 
based systems switchable by electric pulses has grown 
dramatically [5|-|7| leading to important breakthroughs, 
despite a still poor understanding of the switching mech- 
anisms and, in particular, of the joint heat-electron-ion 
transport at the nanoscale. 

Molecular dynamics simulations [8|, |9| of interacting 
oxygen vacancies have demonstrated the formation of va- 
cancy filaments and clusters characterized by a certain 
degree of orientational order. This picture is in qualita- 
tive agreement with recent measurements [10| , which hint 
at the existence of localized conducting channels. More- 
over, the formation of nanoswitches (buried atomic junc- 
tions) between large clusters of vacancies would be con- 
sistent with the observed quantum conductance and the 
unusual resistance-noise spectrum in TaO^; memristors 
[9| . The question then arises to what extent nanoswitches 
formed in oxide memristors are responsible for resistive 
switching and how to effectively control switching in such 
nano j unctions . 

In this Letter, we consider the simplest possible model 
of joint electron- ion- heat dynamics for a nanoswitch with 
electrons tunneling through a bridging nanoparticle (ei- 
ther an ion or a vacancy) subject to thermal fluctuations. 
Electric current through the switch produces Joule heat- 
ing, which results in local temperature relaxation, de- 
scribed by Newton's cooling law. Surprisingly, we demon- 
strate that this system admits two distinct resistive states 




FIG. 2: Slow temperature relaxation: (a) Temperature de- 
pendence of the memristor resistance, R{T) — {R{t,T)} sim- 
ulated by using Eq. ([ij at fixed T. Black dots correspond 
to the stationary points f{Ts) — nTs/qo- Inset: blow-up of 
R{T) at low T. (b) Simulations of coupled Eqs. JU and ^ 
for a current driven memristor and a symmetric harmonic po- 
tential U{x) [Fig. [ijb)] with loq = 2. The hysteretic loop of 
f vs. [qoY oc I occurs only for non-zero currents /. (c) 
Upon shifting the minimum of the parabolic potential, xq 
off-center [Fig[llc)], hysteresis is gradually suppressed and 
negative differential resistance sets on; xq is increased from 
0.3 to 0.9 by steps of 0.1; (d) Resistance hysteretic loops 
for a switch with asymmetric two-well potential, see text, 
with (a, /3, 7) = (4, 11.6, 8) (blue curve) and (4, 18.8, 16) (red 
curves). Other simulation parameters are: time integration 
step At = 0.005, averaging time tmax = 5 x 10*, A = 0.1, 
cjo = 2, and k — 0.01. 



even if the nanoparticle has only one mechanical attractor 
(e.g., the particle moves in a potential with a single min- 
imum). Moreover, if the particle is confined by a poten- 
tial with, say, two minima, these two resistive states can 
survive even at zero applied current (or voltage), so that 
the junction could serve as a non- volatile resistive switch. 
Furthermore, by applying current-voltage pulses such a 
memristor can be made switch into the desired resistive 
state, almost deterministically, although the particle is 
subject to random forces, only. Therefore, the resistance 
of this switch can be cycled through its hysteresis loop 
in apparent analogy with the elementary mechanism of 
magnetic memory devices. In this regard, our model ex- 
hibits all key features of a resistive memory device. 

Model. — Let us consider a nanoparticle diffusing be- 
tween two electrodes subject to thermal noise (Fig. [TJi). 
An electric current through the switch sets on when the 
electrons tunnel from one electrode to the other via the 
nanoparticle. Here, we assume electron relaxation to be 
much faster than heat diffusion and the nanoparticle dy- 
namics, thus, ignoring any shuttling effects [llj. There- 
fore, the instantaneous resistance of the device, R{x), is 
a function of the particle location, x, whereas its tem- 



perature, T{x), results from the balance between Joule 
heating, which depends on R{x), and temperature re- 
laxation. This simplified electron-particle-heat dynamics 
can be formulated as 



dx/dt + dU/dx = V2T^, 
dT/dt = q{R{x)) - kT, 



kl<l, 



(1) 
(2) 



where x is the dimensionless coordinate of the nanopar- 
ticle, which is confined by the binding potential U{x) 
between the two electrodes at x = ±1 (normalization 
condition). The random force S,{t) represents a Gaus- 
sian, zero-mean, white noise with (C(O)^(i)) = S{t). In 
Eq. (m the time dependence of the temperature is de- 
scribed by the standard Newton's cooling law with Joule 
heat source Q{R) = Cq{R), C and k being the junction 
heat capacitance and the temperature relaxation rate, 
respectively. The tunneling resistance between the first 
electrode and the particle and the particle and the sec- 
ond electrode are, respectively, Ri ex exp[(l — x)/X\ and 
R2 (X exp[(a; + 1)/A], where A is the normalized [1^ tun- 
neling length. The total switch resistance can be rewrit- 
ten as R = Ri + R2 = -Ro'': where r = cosh(a;/A) is a 
dimensionless resistance and Rq is the junction resistance 
when the particle is located at the midpoint between the 
electrodes, a; = 0. If the particle is driven by an electrical 
current /, Joule's law states that Q = PR, while in the 
case of voltage driven switches, Q ~ V"^ /R. Correspond- 
ingly, in compact notation, q — qor^^ , with n = ±1 and 
go = I^Ro/C or go = V^/RqC, respectively, for current 
or voltage driven switches. The most interesting working 
regimes occur when the temperature relaxation is either 
slow, K <C 1, or fast, k S> 1. 

Slow temperature relaxation. — For a qualitative anal- 
ysis in the slow relaxation regime, the resistance can be 
averaged with respect to the fast thermal fluctuations of 
the particle, {R{x)) = RiT) = Rof{T), so that Eq. 1^ 
can be reduced to 



dT/dt = qQ[f{T)]" -kT, 



(3) 



with stationary points T — Tg satisfying the identity 
Ts = qf)[f{Ts)f / K. 

We start our analysis by assuming that the nanopar- 
ticle moves in a simple harmonic potential U ~ WqX^/2 
[Fig. [IJb)]. For a voltage driven memristor {n = —1), 
there is only one stationary point T^ and no hysteresis 
is observed. In contrast, for a current driven memris- 
tor (n = 1), there are either one (stable) or three (two 
stable and one unstable) stationary points. The equa- 
tion nTs/qo = ^[Ts) is graphically solved in Fig. [D^a). 
By simulating Eq. ([1]) at constant T, we determined 
the curves f{T) (red solid curve). For ^q < go < ^o 
[(7g and gp correspond to the upper and lower dashed 
straight lines in Fig. HJa)], a straight line KT/qo (dot- 
ted) crosses f{T) in three points; for all other go the 



straight line and the curve f{T) cross in one point, only. 
The profile of the curve f (T) can be sketched analytically 
by noticing that the average {r{x)) must be taken with 
respect to the Boltzmann distribution density p{x) = 
^-ugx /2kT ^ j^^^^^-ujgx jikT f^j. |^| ^ ^ and zero over- 
wise. It follows immediately that the curves f(T) grow 
exponentially like exp[T/(2a;oA^)] at low temperatures, 
T/wq ^ 1, and saturate at Asinh(l/A) for high temper- 
atures, T 3> cJq, in agreement with our simulations [Fig 



[2ja)]. Therefore, within the interval %' < qo < Qq , one 
can expect two stationary solutions of Eqs. ([T]) and ([2]) 
and hysteresis of both temperature and resistance, when 
the applied electrical current is cycled. These conclu- 
sions are confirmed by the numerical integration of the 
coupled Eqs. dH) and © for k = 0.01 [Fig. ^h)]. The 
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characteristics curve displayed there, r{qQ ) (note that 

1/2 

^Q oc /), exhibits hysteretic behavior, thus suggesting 
that this system can serve as a memory bit with two logic 
states corresponding to two different resistances. How- 
ever, such a bit is not a non- volatile memory element, 
yet, since the two states exist only for a nonzero current 
(i.e., qo > 0). 

The resistance hysteretic loop gradually shrinks as 
the position xq of the potential minimum of U{x) ~ 
u!q{x — a;o)^/2 shifts away from the sample center, x — 
where r{x) has its minimum. More interestingly, negative 
differential resistance (NDR) was observed for a;o > A. 
Indeed, at relatively large xq [Fig. [ijc)], the resistance 
first decreases with the current [Fig. [D^c), top red curve] 
and then increases for higher currents. With decreasing 
xq, the NDR region shrinks and finally disappears; on 
further lowering xq, resistance hysteretic loops open up 
and expand (Fig. 2(c), two bottom curves). NDR is due 
to the fact that on increasing T the sharp peak of p(x) 
(which is centered around the potential minimum but 
far away from the r(x) minimum) broadens; correspond- 
ingly, on increasing the current the particle is allowed to 
dwell longer in the vicinity of the resistance minimum at 
x = 0. 

In order to use our device as a non-volatile memory 
element, the nanoparticle must be confined to a potential 
U{x) with two wells at least. Below we analyzed in detail 
the case of the potential U = ax'^ /2 + /3x^/3 + '-fx'^/A with 
two minima. Parameters a,/3,"f control locations and 
depth of the minima. Hereafter, we will consider either 
even potentials with /3 = [Fig. [ijb)] or asymmetric 
potentials [Fig. [Ijc)] with left minimum in the middle, 
X = 0, and right minimum at x sa 0.9, i.e., near the r.h.s. 
electrode. 

For slow temperature relaxation, k <C 1, Fig. ^d) 
clearly shows how the resistance hysteresis gets sup- 
pressed with increasing the depth of the potential minima 
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at X « 0.9: The looped blue curves of r vs. ^q evolve 
into univalued curves - a "virgin" branch showing up 
only if the nanoparticle was initially located in the ecu- 
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FIG. 3: Fast temperature relaxation, (a) Spatial probabil- 
ity distribution of the particle in the current driven switch of 
Eqs. fl]), ^ with double-well potential U{x) with (a,/3,7) = 
(1,0,4). The two-peaked (red) and one-peaked curve (blue) 
represent the low (^q — 0.1) and high current (go — 
0.6) regime, respectively (see text); (b) Spectra of the x{t)- 
fluctuations at low and high current as in (a). The power 
laws 1// and 1//^ are drawn as appropriate for reader's con- 
venience, (c) As in (b) but for the r{t) fluctuations. Here, 
the 1// branch of the spectrum is more pronounced due to 
strong nonhnearity of the r{x) function; (d) Resistance cy- 
cling between low and high resistive state for the double-well 
potential [Fig. [He)], with (a,/3,7) = (4,13.4,8). The volt- 
age, [go]^ oc V, was increased from to 3 on the r.h.s. of 
the loop and the current, [go]^ oc I, from to 0.02 on the 
l.h.s. Time-dependence of the particle location (e), tempera- 
ture (f) and resistance (g) when the voltage pulse, q^ = 2.5, 
applied for < f < 1000, is followed by the current pulse, 
[qo]^^'^ = 0.01, for 1000 < t < 2000 (see text). 



tral well. Therefore, current cycling (or current pulse) 
allows switching the system from the low (nanoparticle 
sitting in the central well) to the high resistive state (par- 
ticle being located near an electrode) , but not vice versa! 
Fast temperature relaxation. — Let us consider now the 
case of fast temperature relaxation, k 3> 1. For a qual- 
itative analysis, we can assume that the temperature is 
always stationary with T(x) = qo[r(x)]"' / k. Substituting 
this expression into Eq. ([T]) for the particle diffusion, we 
derive the Stratonovich stochastic equation 

dx/dt + dU/dx = [cosh(x)]"/2[2go/K]^/^C, l^^l < 1- (4) 

This equation can be regarded as a model for a Brow- 
nian particle diffusing through a medium with spatially 
varying temperature (or diffusion coefficient). For a cur- 
rent driven memristor, n = 1, the temperature reaches 
its maximum when the particle is located near the elec- 
trodes, X = ±1, while for the a voltage driven memris- 
tor the highest temperature occurs when the particle sits 



around a: = 0. Therefore, by applying current or voltage 
pulses, we can "heat up " the center of the device or the 
regions near its electrodes, respectively. This effect can 
be invoked to reliably switch the device between low and 
high resistive state, with little regard for the parameters 
of the current or voltage pulses. Thus, a key problem 
in memristor design can be solved based on the present 
proposal. 

We now propose a switching protocol for our model 

memristor with asymmetric potential U having minima 

in the center, x = 0, and close to the right electrode, x ~ 

0.9, namely a = 4, /3 = 13.4, and 7 = 10. Let us assume 

that the system is initially in the low resistive state. With 

1/2 
increasing voltage, Qq oc V, the memristor switches to 

the high resistive state as shown in our simulations [Fig. 

Eljd); right side of the loop driven by a voltage source]. 

Upon decreasing the voltage to zero, the system sets in 

its high resistive state, indicating that the particle is now 

trapped near the right electrode. If we next increase 

1/2 
the current, q^ oc /, the system switches to the low 

resistive state [Fig. Eld); left side of the loop driven by 

current source]. By finally turning the current off, the 

resistance remains set to its low value. The loop being 

almost deterministic, the proposed switching protocol is 

clearly repeatable. 

In order to illustrate the underlying switching dynam- 
ics, we computed x(t), T{t), and R{t) by numerically 
integrating the coupled Eqs. (P) and © and switch- 
ing n from —1 (voltage driven) to 1 (current driven) at 
t = 1000. This mimics the sequence of voltage pulse, 
with duration < i < 1000, followed by a current pulse 
for 1000 < t < 2000 [Fig. [3i;e)-(f)]. Under the action 
of a constant voltage, < i < 1000, the initial low re- 
sistive state results in strong Joule heating, V^/R, and, 
therefore, high temperature [Fig. [3Kf)] and strong fluc- 
tuations of the nanoparticle position [Fig. [31(e)] just af- 
ter applying a voltage pulse. However, as soon as the 
particle hops over the barrier and gets trapped in the 
well near the right electrode, the resistance increases and 
the heat production drops so fast that the particle tem- 
perature and its spatial fluctuations get quenched [Fig. 
[3l(e)-(f) for t < 1000]. No further dynamics was observed 
over time until after the current pulse was applied for 
t > 1000. Therefore, this switching mechanism is robust 
with respect to voltage pulse duration: Indeed, while the 
memristor is driven by the voltage pulse, the particle can- 
not hop back to the minimum at a; = 0, being the tem- 
perature felt by the particle near the electrode very low. 
When a current pulse is applied, 1000 < t < 2000, the 
heat production starts anew because the particle is now 
in a high resistive state and Joule's heating is propor- 
tional to the switch resistance, Q = I^R. As a result, 
the temperature sharply increases and the particle starts 
to fluctuate around the minimum near the electrode [Fig. 
|3Kc) for t just after t = 1000]. Eventually, the particle 



jumps into the central potential minimum where the re- 
sistance is the lowest. At that very moment the heat 
production drops and the particle gets stuck at a; = 0; no 
more changes are expected as the particle temperature is 
now very low. A strong enough current pulse, irrespec- 
tive of its duration, suffices to switch the system back to 
its low resistive state and keep it there. 

Pulse driven switching in our model turns out to be so 
effective because resistance and heat production strongly 
correlate with the particle location and the particle is al- 
ways attracted toward the region with the lowest heat (or 
entropy) production, where it gets trapped. This conclu- 
sion can be easily extended to the case of any memristor 
with two distinct resistive states with resistances i?i and 
i?2 such that i?i <C i?2- As a non-equilibrium thermo- 
dynamic system tends to attain a state of lowest entropy 
production, a voltage pulse will switch the memristor to 
the high resistive state i?2, while a current pulse will 
switch it back to the low resistive state Ri. Thus, our 
conclusions about optimal switching sequence are very 
general and can be readily extended to more elaborated 
memristor models. 

Finally, we stress that the memristor switching mech- 
anism proposed here can be regarded as an instance of 
noise- induced phase transition [13| . Such a stochastic dy- 
namical phenomenon takes place when the minimum of 
the temperature T{x) does not coincide with the mini- 
mum of the binding potential, U{x): The Brownian par- 
ticle tends to dwell about the potential minima at low 
temperatures, and about the temperature minima at high 
temperatures, thus switching location on increasing the 
noise. In order to detect a noise-induced transition in 
our memristor model, we simulated the coupled Eqs. ([T]) 
and ([2]) for a symmetric double well potential U{x) with 
a = l,/3 = 0, and 7 = 4, in the current driven regime 
with K = 25 (fast temperature relaxation) [Fig. |3lja)]. 
For weak currents, the stationary temperature T(x) is 
low and the particle is mostly localized around the po- 
tential minima; the corresponding x distribution, p{x), is 
doubly peaked (red curve). At higher currents, the tem- 
perature increases until the particle is pushed out of the 
potential wells toward the center of the sample, where 
the potential has a maximum, but the temperature is 
much lower. In the fast relaxation regime of Eq. (|4]) such 
a transition is well described by stationary probability 
distribution p{x) = [J\f/,yT{x)] exp[- J^ dxU'{x)/T{x)], 
where U' = dU/dx and A/" is a normalization constant. A 
current increase is accompanied by a temperature surge, 
so that the minimum of the distribution at a; = can 
change to a maximum [l3| . More interestingly, the spec- 
tra of the resistance and position fluctuations (also called 
noise) show a qualitative change in correspondence with 
such a transition: The standard 1//^ spectra break up 
into a low-frequency 1// branch and a high-frequency 
1//^ tail [Fig. [2Ib),(c)]. The unusual experimental noise 
spectra reported in Ref. [9| can thus be interpreted as 



the signature of a noise-induced transition occurring in 
TaOa; memristors. 

Conclusion. — The resistive switch model discussed 
here is expected to ensure the much sought-for proper- 
ties of controhabihty and rehabihty one needs to design 
and prototype switchable nanomemristors for future flash 
memory devices. On the other hand, based on numerical 
simulation, we claim that its response to external cur- 
rent and voltage signals reproduces well a wide variety of 
experimental observations reported in the literature for 
real memristors. Moreover, the predicted phenomenon of 
negative differential resistance can considerably enhance 
the functionality of these nanostructurcs. 

Wc acknowledge many fruitful discussions with R.S. 
Williams, who strongly supported this work. S.E.S. 
also acknowledges The Leverhulme Trust for partial sup- 
port of this research; P.M. aknowledges partial sup- 
port from the European Commission, grant No. 256959 
(NanoPower). 



[1] L. Venema, Nature 479 309 (2011). 

[2] B. Dipert, M. Levy Designing with Flash Memory (Anna- 
books, San Diego, 1993). 

[3] J. Palmer, Science and technology re- 

porter, BBC News, 16 May 2011: 

http://www.bbc.co.uk/news/science-environment-13392857 

[4] G. Dearnaley, A.M. Stoneham, and D.V. Morgan, Rep. 
Prog. Phys. 33, 1129 (1970). 

[5] D.B. Strukov, G.S. Snider, D.R. Stewart, and R.S. 



Williams, Nature 453, 80 (2008); K. Szot, W. Speier, 
G. Bihlmayer, and R. Waser, Nature Materials 5, 312 
(2006). 
[6] J.J. Yang et al, Nat. Nanotechnol. 3, 429 (2008) 
[7] J.J. Yang, F. Miao, M.D. Pickett, D.A.A Ohlberg, D.R. 
Stewart, C.N. Lau, and R.S. Williams, Nanotechnology 
20, 215201 (2009). 
[8] S.E. Savel'ev, A.S. Alexandrov, A.M. Bratkovsky, and 
R.S. Williams, Nanotechnology 22, 254011 (2011); S.E. 
Savel'ev, A.S. Alexandrov, A.M. Bratkovsky, and R.S. 
Williams, Appl. Phys. A 102, 891 (2011); S.E. Savel'ev, 
A. S. Alexandrov, A.M. Bratkovsky, and R.S. Williams, 
Appl. Phys. Lett. 99, 053108 (2011). 
[9] Wei Yi, S.E. Savel'ev, G. Medeiros-Ribeiro, J.J.Yang, 
A.M. Bratkovsky, and R.S. Williams, to be published. 

[10] J.J. Yang et al., Nat. Nanotechnol. 3, 429 (2008); 
F. Miao, J.J. Yang, J. P. Strachan, D. Stewart, R.S. 
Williams, and C.N. Lau, Appl. Phys. Lett. 95, 113503 
(2009); R. Mustermann et al, Phys. Status Sohdi 4, 
16 (2009); J. P. Strachan, M.D. Pickett, J.J. Yang, S. 
Aloni, A.L.D. Kilcoyne, G. Medeiros-Ribeiro, and R.S. 
WiUiams, Adv. Mater. 22, 3573 (2010); D-H. Kwon et 
al, Nat. Nanotechnol. 5, 148 (2010). 

[11] L. Y. Gorelik et al, Phys. Rev. Lett. 80, 4526 (1998); 
N.A. Zimbovskaya, M.R. Pederson, Physics Reports 509, 
1 (2011); A.V. Moskalenko et al, Phys. Rev. B 79, 
241403(R) (2009). 

[12] i.e. we rescaled A by half the distance between the elec- 
trodes. 

[13] W. Horsthemke and R. Lefever, Noise-Induced Transi- 
tions (Springer, Berlin, 1984); C. Festa, L. Fronzoni, 
F. Marchesoni, and P. Grigolini, Phys. Lett. A 102, 95 
(1984). 



